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PREFACE 


This  report  presents  the  results  of  modeling  the  effects  of  electron  density  irregulari¬ 
ties  on  transionospheric  propagation.  The  objective  of  this  effort  has  been  to  provide  a  tool 
for  determining,  in  real  time,  the  expected  effects  of  such  irregularities  on  the  performance 
of  a  given  transionospheric  monitoring  system.  While  these  effects  can  take  several  forms, 
we  concentrate  here  on  the  variation  in  the  phase  of  the  signal  due  to  scattering  caused  by 
such  irregularities. 

In  this  report,  we  shall  be  mainly  interested  in  remote  locating  systems.  Such 
systems  rely  on  accurate  measurements  of  the  phase  difference  between  separated  receive 
antennas  to  provide  estimates  of  apparent  angular  location  (azimuth/ elevation).  Fluctua¬ 
tions  in  this  phase  difference  due  to  scattering  caused  by  electron  density  irregularities 
result  in  a  consequent  fluctuation  in  the  apparent  location  of  the  transmitter.  System 
performance  can  be  seriously  degraded  if  the  correlation  in  the  phase  fluctuation  between 
antennas  is  low  or  the  mean  square  fluctuation  in  phase  is  large.  In  any  case,  the  accurate 
determination  of  the  phase  fluctuations  provides  an  estimate  of  the  error  for  a  given  loca¬ 
tion  estimate  and  thus  provides  important  information  for  the  system  user. 

Another  feature  that  must  be  addressed  under  strong  multiple-scattering  conditions 
is  the  reduction  of  correlation  in  the  intensity  of  the  signal.  Under  such  scattering  condi¬ 
tions,  defocusing  by  large  scale*  irregularities  can  cause  deep  small-period  fades  in  signal 
strength.  For  sufficiently  strong  scattering,  these  signal  fades  may  have  a  spatial  extent  on 
the  order  of  the  size  of  the  locating  system  (tens  of  meters).  This  loss  of  signal  makes  the 
system  unreliable  at  best  and,  in  some  cases,  inoperable. 

This  area  of  investigation  is  not  completely  new;  there  exists  a  large  body  of  related 
work.  What  is  different  about  this  effort  is  the  extension  to  low  VHF  and  HF  operational 
frequencies.  Previous  work  has,  for  the  most  part,  dealt  with  higher  frequencies  where,  in 
general,  the  relatively  simple  single-scatter  models  are  applicable.  Also,  at  these  frequen¬ 
cies,  one  can  generally  ignore  the  refractive  effects  of  the  background  ionization.  This 
simplifies  the  implementation  of  the  models  considerably. 

At  lower  frequencies,  strong  multiple  refractive  scattering  by  large  scale  irregulari¬ 
ties  becomes  the  rule  rather  than  the  exception.  Also,  refractive  effects  of  the  background 
ionosphere  become  extremely  important  and  must  be  included  for  accurate  modeling.  Thus 
applicable  models  become  much  more  complicated  and  their  use  in  a  real-time  situation 
doubtful.  However  over  the  last  decade  or  so.  Professor  H.G.  Booker  of  the  University  of 
California  at  San  Diego  has  developed  a  model  suitable  in  these  situations.  This  model  is 
described  and  used  in  this  report  to  provide  a  means  of  calculating  the  effects  described 
above. 


*In  this  report  the  term  “scale”  will  always  refer  to  wavelength  (A)  divided  by  2  n\ 
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INTRODUCTION 


In  this  report,  we  consider  the  propagation  of  electromagnetic  waves  within  a 
bounded  plasma  medium.  If  the  characteristics  of  such  waves  are  known  before  their  pene¬ 
tration  into  such  a  region,  measuring  the  same  characteristics  after  the  waves  exit  is  a  direct 
means  of  determining  the  properties  of  the  plasma  medium  itself.  For  example,  the  devia¬ 
tion  of  ray  direction  for  a  pulse  of  suitable  carrier  frequency  after  traversing  the  earth’s 
ionosphere  provides  a  means  of  determining  the  magnitude  of  electron  density  gradients 
within  the  ionosphere. 

Another  example  is  the  modulation  of  the  wavefront  which  often  occurs  when 
a  wave  traverses  the  ionosphere.  This  modulation  is  caused  by  irregularities  of  electron 
density  which  exist  within  the  ionosphere.  The  variation  in  phase  across  the  wavefront 
produces  interference,  and  eventually,  amplitude  fluctuations  at  a  distant  receiver.  One  way 
of  determing  the  amplitude  and  scale  of  the  irregularities  within  the  ionosphere  is  to 
measure  the  magnitude  and  scale  of  these  amplitude  fluctuations.  This  is  the  phenomenon 
known  as  amplitude  scintillation  and,  along  with  the  phase  fluctuations,  known  as  phase 
scintillation,  is  the  subject  of  this  report. 

For  transionospheric  propagation,  the  effects  of  the  scattering  can  take  several 
forms,  depending  on  the  characteristic  of  the  wave  being  measured.  For  transionospheric 
communication  systems,  the  effects  of  scattering  can  be  significant.  The  deep  high- 
frequency  fading  of  signal  levels  having  a  duration  of  minutes  or  even  hours  can  cause 
extended  periods  for  which  communication  may  be  impossible. 

Scattering  can  also  affect  burst  communication  systems,  which  require  coherent 
signals  over  a  large  bandwidth.  Frequency  dependence  of  the  scattering  causes  significant 
dispersive  effects,  which  can  drastically  reduce  effective  signal  bandwidth. 

Locator  systems,  which  use  multiple  receiving  antennas  for  remote  detection,  are 
also  greatly  affected.  Such  systems  rely  on  differences  in  the  phase  of  received  signals  at 
each  antenna  to  determine  the  location  of  the  transmitter.  Fluctuations  in  the  phase  of  the 
signal  caused  by  the  scattering  process  cause  apparent  variations  in  the  location  of  the 
transmitter.  Under  very  strong  scattering  conditions,  the  phase  fluctuations  can  become 
large  enough  to  render  such  locator  systems  incapable  of  operation. 

Several  methods  are  available  for  modeling  the  effects  of  scattering  on  transiono¬ 
spheric  propagation.  Methods  employed  differ,  depending  on  whether  the  scattering  is 
weak  or  strong.  If  the  scattering  is  weak,  which  we  will  define  more  rigorously  later,  the 
weak  or  single-scatter  modeling  methods  are  applicable.  These  include  the  Born  Approxi¬ 
mation  (Ishimaru,  1978),  the  Rytov  or  smooth  perturbation  method  (Tatarskii,  1971; 
Ishimaru,  1978),  and  the  phase-changing  screen  method  of  Booker  (Booker  et  al .,  1950). 
These  methods  are  usually  adequate  at  very  high  frequencies  or  under  extremely  quiescent 
ionospheric  conditions. 

At  lower  frequencies  oi  under  highly  disturbed  conditions,  strong  or  multiple 
scattering  becomes  the  rule.  Theories  which  describe  the  scattering  under  these  conditions 
are  generally  extremely  complicated  and  are  therefore  not  amenable  for  use  in  an  opera¬ 
tional  environment.  Two  of  the  currently  popular  methods  of  this  type  in  use  are  based  on 
the  parabolic  approximation  to  the  wave  equation.  These  are  the  path  integral  method 
(Dashen,  1979)  and  the  moment  method  (Ishimaru,  1978). 

Another  much  simpler  method  is  the  multiple  refractive  scatter  model  developed 
by  Booker  (Booker,  MajidiAhi,  1981).  In  this  method,  the  multiple-scattering  region  is 
replaced  by  a  single  “deeply  modulated”  phase-changing  screen.  This  amounts  to  replacing 
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a  multiple  small-angle  scatter  situation  with  a  single  “large”  angle  one.  At  first  glance,  the 
validity  of  such  an  approximation  is  questionable.  However,  it  has  been  shown  (Booker 
el  al.,  1985)  that  as  long  as  the  centrally  located  phase  screen  (1)  has  effective  irregularity 
sizes  perpendicular  to  an  incident  ray  which  are  identical  to  those  in  the  medium,  and  (2) 
produces  the  same  total  mean-square  value  of  phase  fluctuations  as  the  actual  medium,  the 
theory  produces  results  substantially  the  same  as  the  more  elaborate  theories  mentioned 
above.  The  relative  simplicity  of  the  multiple  refractive  scatter  model  also  makes  it  attrac¬ 
tive  for  use  in  an  operational  environment.  For  these  reasons,  this  model  was  chosen  for 
this  project  and  is  the  model  described  in  this  report. 

Any  model  used  to  predict  the  effects  of  F-region  electron  density  irregularities 
on  transionospheric  propagation  requires  an  accurate  description  of  the  properties  of  the 
irregularities  themselves.  Fortunately,  a  large  literature  exists  detailing  the  results  of  exper¬ 
iments  designed  to  measure  these  properties.  In  the  next  section,  we  give  a  general  overview 
of  the  findings. 

F-REGION  ELECTRON  DENSITY  IRREGULARITIES 

A  great  deal  of  research  has  been  done  on  the  mechanisms  which  produce  electron 
density  irregularities  in  the  ionosphere.  Generally  speaking,  it  has  been  found  that  the 
experimentally  determined  properties  of  the  irregularities  are  consistent  with  various 
plasma  wave  instability  processes.  These  processes  differ  on  a  global  scale,  so  that  the 
mechanism  thought  to  produce  equatorial  “bubble”  irregularities,  for  example,  is  different 
from  that  which  is  supposed  to  produce  auroral  “blob”  irregularities.  This  variability  of 
causative  mechanisms  makes  difficult  the  development  of  a  unified  global  picture  of  the 
generation  of  electron  density  irregularities. 

A  discussion  of  the  processes  for  the  generation  of  ionospheric  irregularities  is 
outside  the  scope  of  this  report.  Readers  interested  in  further  information  are  directed  to 
Fejer,  Kelley,  1980;  Heelis,  1987;  and  Tsunoda,  1988.  In  this  report,  we  simply  assume  the 
existence  of  the  irregularities  in  the  F-region  ionosphere  and  attempt  to  describe  their 
general  physical  properties.  Later  we  will  relate  these  properties  to  the  observed  properties 
of  the  scattering  phenomena  they  produce. 

Ground-based  scintillation  measurements  are  one  means  of  probing  the  irregular 
structure  of  the  ionosphere.  Among  other  important  measurement  techniques  are  ground- 
based  polarization  rotation  measurements  (Yeh  et  al.,  1979),  in  situ  measurements  by 
rocket  and  orbiting  satellites  (Singh,  Szuszczewicz,  1984;  McClure,  Hanson,  1973;  Clark, 
Raitt,  1976),  and  radar  backscatter  measurements  (Woodman,  LaHoz,  1976).  All  these 
techniques  have  been  used  to  produce  a  general  description  of  the  physical  properties  of 
F-region  electron  density  irregularities  for  scale  sizes  from  several  meters  to  hundreds  of 
kilometers. 

One  important  characteristic  of  ionospheric  plasma,  which  we  will  not  discuss  in 
depth  in  this  report,  is  its  motion  with  respect  to  the  earth.  This  motion  is  caused  by  forces 
produced  by  various  natural  mechanisms,  among  which  are  neutral  atmospheric  winds, 
neutral  atmospheric  acoustic  gravity  waves,  and  the  interaction  of  the  earth’s  magnetic  field 
with  induced  electric  fields.  These  forces  cause  the  plasma,  and  hence  the  irregularities,  to 
drift  in  complicated  patterns. 

In  all  experimental  systems  designed  to  measure  the  spatial  properties  of  electron 
density  irregularities,  the  data  are  actually  determined  in  the  form  of  a  time  series.  For 
example,  a  single  ground-based  receiver  in  a  transionospheric  scintillation  experiment 
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measures  the  changes  in  the  signal  as  a  function  of  time.  This  temporal  information  is 
multiplied  by  an  “effective”  velocity  of  the  irregularities  to  provide  information  about  their 
spatial  characteristics.  This  effective  velocity  is  determined,  in  any  given  experiment,  by  any 
or  all  of  the  natural  mechanisms  listed  above.  Thus,  it  is  possible  to  describe  the  properties 
of  electron  density  irregularities  in  either  temporal  or  spatial  terms.  In  this  report,  we 
discuss  the  irregularities  in  spatial  terms,  as  is  customary  in  the  literature. 

The  properties  of  the  irregularities  important  in  determining  the  nature  of  the 
scattering  are  their  amplitude  and  their  spatial  power  spectral  density  or  power  spectrum. 
The  amplitude  determines  the  strength  of  the  interaction  between  the  irregularities  and  an 
electromagnetic  wave.  The  power  spectrum  determines  the  spatial  extent  of  the  interaction. 
We  shall  treat  each  separately  below. 

AMPLITUDE  CHARACTERISTICS 

Most  of  the  current  knowledge  concerning  the  fluctuation  amplitude,  A N/N,  of 
electron  density  irregularities  has  been  gained  through  the  use  of  rockets  and  orbiting  satel¬ 
lites.  These  in  situ  methods  only  produce  information  in  one  spatial  dimension,  the  direc¬ 
tion  of  travel  of  the  platform,  so  a  complete  three-dimensional  description  requires  the 
combination  of  results  from  various  experiments. 

Generally,  the  irregularity  amplitudes  that  are  observed  differ  between  high  lati¬ 
tudes,  mid-latitudes,  and  equatorial  latitudes.  At  mid-latitudes,  the  F-region  ionosphere 
seems  to  be  remarkably  smooth.  An  analysis  of  data  from  the  satellite  Ogo  6  (perigee  ~400 
km,  apogee  ~1 100  km,  inclination  ~82  degrees)  by  McClure  and  Hanson  (1973)  shows 
average  root-mean-square  (rms)  values  of  amplitude  of  about  0.1%  for  irregularities  with 
wavelengths  less  than  100  km.  (Irregularity  sizes  are  usually  given  in  either  wavelength  or 
scale.)  This  value  is  also  observed  in  nighttime  when  scintillation  activity  is  usually  most 
intense  (McClure,  Hanson,  1973). 

A  more  extensive  data  set  collected  by  the  satellite  ESRO-4  (perigee  ~240  km, 
apogee  ~  1 1 77  km,  inclination  ~91  degrees)  was  analyzed  by  Clark  and  Raitt  (1976).  These 
authors  organized  approximately  1.5  years  of  rms  amplitude  data  into  contour  plots  of 
location  vs.  local  time.  Due  to  the  precession  of  the  orbital  plane,  it  was  possible  to  order 
their  data  according  to  altitude.  The  results,  valid  for  irregularity  wavelengths  between 
about  7  m  and  7  km,  show  that  for  magnetically  quiet  periods  (Kp  <  3),  average  values  of 
rms  amplitude  are  less  than  or  about  2%.  This  is  for  amplitude  values  averaged  over  a 
height  range  of  400-1000  km.  Again,  this  value  is  for  both  day  and  night  periods  (Clark 
and  Raitt,  1976). 

In  the  restricted  height  range  of  400-600  km  (i.e. ,  just  above  the  F-region  peak) 

®  subauroral  zone  structure  tends  to  merge  with  mid-latitude  structure  at  nighttime  to 

produce  average  rms  amplitudes  as  high  as  4%.  This  structure  also  extends  to  the  equator¬ 
ial  region  where  contours  from  the  northern  and  southern  hemispheres  join.  This  produces 
the  only  evidence  of  significant  equatorial  irregularities  (~4%)  seen  by  Clark  and  Raitt 
(1976).  Similar  results  were  found  for  magnetically  disturbed  periods  although,  in  this  case, 
nighttime  subauroral  structure  is  evident  (~2%)  at  mid-  and  equatorial  latitudes  even  when 
averaged  over  the  full  height  range  of  400-1000  km  (Clark.  Raitt,  1976). 

One  phenomenon  observed  at  all  latitudes,  but  especially  in  the  equatorial  region, 
is  “spread-F.”  The  term  “spread-F"  is  used  to  describe  a  scattering  phenomenon  in  which 
F-region  ionograms  show  echoes  which  are  spread  over  a  range  of  heights  or  frequencies. 
When  a  range  of  frequencies  near  the  peak  frequency  shows  returns  from  many  height 
ranges,  the  condition  is  often  termed  “range  spread-F”  and  is  associated  with  strong  scintil¬ 
lation  activity  (Aarons,  1982). 
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The  spread-F  phenomenon  has  been  the  subject  of  much  research  since  it  was  First 
observed  in  the  mid-30s  (Booker,  Wells,  1934).  The  equatorial  occurrence  of  the  phenom¬ 
enon  is  thought  to  be  caused  by  large  scale  ( —  1 00  km)  bubble  irregularities  observed  in 
postsunset  hours.  These  regions  of  depleted  electron  density  are  believed  to  be  generated  by 
bottomside  instability  processes,  where  they  produce  the  spread-F  signature  on  ionograms 
(Fejer,  Kelley,  1980).  Once  formed,  the  bubbles  rise  to  heights  above  the  F-layer  peak, 
where  they  produce  transionospheric  scintillation  effects  (Fejer,  Kelley,  1980).  The  scintil¬ 
lation  is  believed  to  be  caused  by  smaller  scale,  large-amplitude  (~100%)  irregularities 
formed  at  the  steep  density  gradients  at  the  edges  of  the  bubble.  This  phenomenon,  while 
somewhat  localized  in  time  and  space,  produces  the  most  intense  scintillation  observed 
(Aarons,  1982). 

At  high  latitudes,  the  location  of  the  auroral  oval  provides  a  scintillation  boundary 
between  mid-latitude  and  auroral  zone  irregularity  structure.  The  analysis  of  Clark  and 
Raitt  (1976)  shows  that  at  all  height  levels  from  400  to  1000  km,  the  equatorial  boundary  of 
the  auroral  scintillation  region  reaches  a  geomagnetic  latitude  of  ~40  degrees  in  the  night¬ 
time  sector.  Within  this  nighttime  subauroral  region,  these  authors  find  that  almost  60%  of 
the  time,  rms  values  of  amplitude  are  greater  than  5%;  and  almost  100%  of  the  time,  they 
are  greater  than  2%.  Peak  values  in  this  region  reach  8%- 10%  at  all  height  ranges.  This 
behavior  extends  poleward  of  the  auroral  zone,  with  only  a  small  decrease  in  observed 
intensity  (Clark,  Raitt,  1976). 

In  summary,  the  experimental  data  show  that  irregularity  amplitudes  differ  accord¬ 
ing  to  location  and  local  time.  At  high  latitudes,  the  irregularities  are  generally  of  largest 
amplitude,  especially  in  the  region  of  the  nighttime  auroral  zone.  These  large-amplitude 
irregularities  continue  only  slightly  diminished  across  the  polar  region. 

At  mid-  and  equatorial  latitudes,  irregularities  are  generally  much  less  intense. 
However,  at  heights  close  to  the  F-region  peak,  nighttime  subauroral  zone  structure  tends 
to  reach  equatorward.  This  produces  the  only  significant  average  irregularity  intensities  in 
this  region. 

In  a  more  localized  sense,  large-amplitude  irregularities  are  often  seen  in  coinci¬ 
dence  with  nighttime  spread-F.  While  spread-F  can  occur  at  any  latitude,  it  is  most  often 
seen  at  equatorial  latitudes,  where  it  produces  the  most  intense  scintillation  activity. 

SPECTRAL  PROPERTIES  OF  IRREGULARITIES 

Electron  density  irregularities  in  the  F-region  ionosphere  have  been  observed  over 
a  large  range  of  scale  sizes,  from  several  meters  to  thousands  of  kilometers.  The  largest  of 
these  irregularities  are  associated  with  tidal  motions  and  planetary  scale  wave  motions 
(Livingston  et  al.,  1981),  while  the  smallest  are  associated  with  the  turbulent  decay  of 
intermediate  scale  structure^  Between  these  extremes  lies  a  scale  range  for  which  the  spatial 
power  spectrum  can  be  characterized  as  power  law.  The  irregularities  responsible  for  the 
scintillation  observed  on  the  earth  lie  largely  within  this  range. 

The  extremes  of  the  power  law  scale  range  are  determined  by  two  particular  scale 
sizes.  At  the  small  scale  (large  wave  number)  end,  the  “inner  scale”  is  thought  to  lie  at 
approximately  the  ion  gyroradius  (~3  m)  (Booker,  Ferguson,  1978;  Woodman.  Basu,  1978). 

At  the  large  scale  (small  wave  number)  end  of  the  spectrum,  there  is  some  debate 
as  to  the  location  of  the  “outer  scale.”  Some  workers  believe  that  no  identifiable  largest 
scale  size  irregularity  exists.  They  argue  that  the  large  scale  irregularities  merge  smoothly 
into  nonirregular  “trends,”  such  as  the  diurnal  effects  of  solar  radiation,  and  so  a  separa¬ 
tion  scale  between  irregular  and  deterministic  effects  is  impossible  to  define.  In  this  case. 
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the  largest  scale  size  of  importance  in  any  operational  situation  is  determined  by  system 
parameters,  such  as  integration  time  or  bandwidth  of  a  receiver  (Rino,  1979). 

However,  others  argue  that  a  definable  effective  outer  scale  does  exist  but  has  not 
been  detected  to  date  because  of  the  det  rending  technique  that  has  been  used  in  analyzing 
measured  data.  The  detrending  procedure  is  commonly  used  to  remove  a  deterministic 
component  from  a  given  data  set.  For  example,  when  phase  variations  of  the  signal  from 
an  orbiting  satellite  are  analyzed,  the  known  effect  of  the  satellite  motion  is  removed  in 
order  to  isolate  the  variation  due  to  scattering.  However,  if  the  detrending  frequency  is 
higher  than  any  frequency  components  contained  in  the  data,  all  the  lower  frequency 
content  is  lost  in  this  procedure.  It  is  the  low-frequency  component  of  the  data  which 
contains  the  effects  of  the  large  scale  irregularities,  and  therefore  removing  this  compo¬ 
nent  makes  it  impossible  to  identify  an  outer  scale  from  phase  scintillation  data. 

Booker  argues  (Booker,  1979)  that  the  turbulence-like  power  law  spectrum  (see 
below)  which  is  observed  results  from  the  cascade  of  energy  from  the  large  outer-scale 
irregularities  to  smaller  scale  sizes.  The  outer-scale  irregularities  are  produced  by  travelling 
ionospheric  disturbances  (TIDs),  which  seem  to  be  present  to  one  degree  or  another  at 
most  times  in  the  F-region  ionosphere.*  TIDs  are  produced  by  the  interaction  of  neutral 
atmospheric  acoustic  gravity  waves  with  the  ionized  component  of  the  atmosphere  (Hines, 
1960).  The  amplitude  of  the  gravity  waves  at  any  height  is  dependent  on  the  scale  height  of 
the  atmosphere  at  that  height  (Hines,  1960).  This  leads  Booker  to  identify  the  outer  scale  of 
the  irregularity  spectrum  with  the  neutral  atmospheric  scale  height  (Booker,  1979;  Booker. 
Tao,  1987). 

In  order  to  resolve  the  question  of  the  existence  of  a  physical  outer  scale  of  turbu¬ 
lence,  an  experiment  must  be  designed  to  determine  the  spectrum  at  low  frequencies.  Until 
then,  the  question  remains  open. 

Between  the  inner  scale  and  outer  scale,  however  it  may  be  defined,  the  power 
spectrum  of  irregularities  is  characterized  by  the  form  k ~p,  where  k  is  the  wavenumber 
(I /scale  size)  and  p  is  the  spectral  index.  In  general,  the  results  from  in  situ  measurements 
(i.e.,  one-dimensional)  show  that  p  ranges  between  1.5  and  2.5,  with  values  close  to  2  being 
most  common.  This  result  is  found  at  both  high  (Phelps,  Sagalyn,  1976)  and  equatorial 
(Livingston  et  al.,  1981)  latitudes  for  irregularities  from  ~ 1 00  m  to  several  hundred 
kilometers. 

Consistent  with  these  results  for  the  power  spectrum,  the  spectral  index  for  the 
amplitude  spectrum  is  approximately  unity  (McClure,  Hanson,  1973).  Thus,  the  average 
amplitude  <AN/  N>  is  approximately  proportional  to  scale  size  (McClure,  Hanson,  1973). 
Further,  the  actual  spatial  gradients  of  electron  density  are  approximately  independent  of 
scale  size,  depending  only  on  the  density  of  the  ambient  plasma. 

Some  more  recent  results  have  indicated  a  systematic  decrease  in  the  spectral 
index  for  strong  perturbation  strengths  (Livingston  et  al .,  1981).  Similar  effects  have  been 
observed  in  phase  scintillation  data  taken  under  strong  scatter  conditions  (Livingston  et  al., 
1981).  Thus,  in  some  modeling  situations,  it  may  be  appropriate  to  consider  two  values  of 
the  spectral  index. 


*A.  Paul  of  NOSC,  private  communication,  1988. 
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The  general  spectral  properties  described  above  have  been  summarized  by  Booker 
(1979)  and  are  shown  in  Fig.  1.  This  figure  actually  represents  the  inferred  two-dimensional 
spectrum  for  irregularities  isotropic  in  planes  perpendicular  to  the  earth’s  magnetic  field. 
For  the  two-dimensional  spectrum,  the  spectral  index  is  p  +  1.  Also  indicated  in  Fig.  1  is 
the  portion  of  the  spectrum  responsible  for  the  various  scattering  phenomena  observed. 

In  situ  data  provide  spectral  information  in  only  one  spatial  direction.  A  full  three- 
dimensional  description  of  the  electron  density  irregularities  requires  multitechnique  exper¬ 
iments.  From  intensity  and  phase  scintillation  measurements  of  transionospheric  signals,  it 
has  been  found  that  the  electron  density  irregularities  are  structured  by  the  earth’s  magnetic 
field.  Results  from  the  Wideband  satellite,  for  example,  confirm  eailier  results  that  the  elec¬ 
tron  density  irregularities  are  aligned  along  the  earth’s  magnetic  field  lines  (Rino,  Living¬ 
ston,  1982).  At  equatorial  latitudes,  the  alignment  is  in  the  form  of  rod-like  structures  elon¬ 
gated  along  the  magnetic  field  lines.  Analysis  of  Wideband  satellite  data  (Secan,  Fremouw, 
1983)  shows  that  for  irregularities  up  to  about  20-km  cross-field  wavelength,  the  ratio  of 
the  scale  size  along  the  field  to  across  the  field  is  about  30:1.  This  ratio  is  reduced  as  the 
latitude  increases  to  about  8:1  (Secan,  Fremouw,  1983). 

At  auroral  latitudes,  the  Wideband  phase  scintillation  data  show  an  increase  in 
the  mean-square  fluctuation  of  phase  when  the  satellite  line-of-sight  lies  within  the  local 
L  shell.  This  is  found  to  occur  even  when  the  satellite  is  not  in  the  same  meridian  as  the 
receiving  site  (Rino  et  al,,  1978).  This  phenomenon  has  been  interpreted  as  implying  an 
elongation  of  the  irregularities  in  the  magnetic  east-west  direction  within  an  L  shell  (Rino 
et  al.,  1978).  The  observed  extent  of  this  elongation  along  a  second  axis  is  a  more  compli¬ 
cated  function  of  local  magnetic  time  and  of  the  position  of  the  line-of-sight  ionospheric 
penetration  point  within  the  auroral  region  (Secan,  Fremouw,  1983). 

The  results  quoted  above  are  based  on  ground-based  measurements  of  the  phase 
scintillation  on  signals  from  the  Wideband  satellite.  These  data  were  routinely  detrended  at 
a  period  of  10  seconds  before  analysis  (Fremouw  et  al.,  1978).  Thus,  the  results  are  strictly 
valid  for  irregularities  with  scale  sizes  of  a  few  kilometers  or  less.  Evidence  indicates  that 
large  scale  irregularities  (~  outer  scale  and  larger)  are  approximately  magnetically  isotropic 
(Booker,  1979). 

In  summary,  the  spectral  properties  of  electron  density  irregularities  are  basically 
independent  of  global  location.  Results  from  different  experimental  methods  all  indicate  a 
power  law  spectrum  with  a  one-dimensional  spectral  index  near  2. 

Irregularities  up  to  about  3  km  are  observed  to  be  structured  by  the  earth’s 
magnetic  field.  This  structure  takes  the  form  of  rod-like  irregularities  at  equatorial  and 
mid-latitudes.  At  auroral  zone  latitudes,  the  rod-like  structure  gives  way  to  a  sheet-like 
structure,  where  elongation  is  within  the  local  L  shell. 

In  the  next  section,  we  describe  the  multiple  refractive  scatter  model.  This  will 
require  an  analysis  of  the  phase-changing-screen  modeling  technique.  First,  however,  we 
provide  an  intuitive  description  of  the  multiple-scattering  process  and  how  it  differs  from 
weak  scattering. 
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Log  (Spectral  Density  of  Irregularities) 


Horizontal  Scale  x 

Scale  (km)  Magnetic  Reid  (m) 

1000  100  10  1  100  10  1  0.1  0.0 1 


Figure  1.  Representation  of  the  power  spectrum  of  iono¬ 
spheric  irregularities.  Portions  of  the  spectrum  responsible 
for  various  scattering  phenomena  are  indicated  on  the 
figure  (after  Booker,  1979). 
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EFFECT  OF  ELECTRON  DENSITY  IRREGULARITIES 
ON  TRANSIONOSPHERIC  PROPAGATION 

Consider  a  plane  electromagnetic  wave  of  frequency / and  wavelength  X  incident 
at  an  arbitrary  angle  on  an  extended  plasma  medium.  If  /is  larger  than  the  peak  plasma 
frequency,  the  wave  will  propagate  into  the  plasma.  Assume  that  this  region  contains  a 
random  distribution  of  electron  density  irregularities  with  sizes  ranging  from  about  kjlrr 
to  some  outer  scale  size  which  is  much  greater  than  k/ 2n.  As  the  wavefront  propagates 
within  the  medium,  the  variation  of  refractive  index  it  encounters  causes  sections  of  the 
wavefront  to  advance  with  respect  to  the  main  wave  and  other  sections  to  recede.  This 
results  in  a  random  “corrugation”  of  the  initially  plane  wavefront  which  we  can  describe  in 
terms  of  the  rms  variation  in  phase,  <A d>2>,/5,  across  the  wavefront.  As  the  wave  continues 
to  propagate,  diffraction  effects  cause  amplitude  variations  to  develop  across  the  wave- 
front.  The  distance  the  wave  must  travel  before  amplitude  variations  develop  is  dependent 
on  <A d>2>'/:,  X,  and  the  scale  size  of  the  irregularities. 

Let  A L  be  the  depth  of  penetration  of  the  wavefront  within  the  medium.  When  the 
wavefront  first  enters  the  medium,  the  rms  variation  of  phase  across  the  wavefront  is  very 
small,  and  hence  any  resulting  amplitude  variations  are  also  very  small.  For  the  ionosphere, 
the  density  irregularities  one  normally  encounters  have  a  small  enough  amplitude  that  even 
for  A L  on  the  order  of  the  outer  scale  size  of  irregularities,  the  rms  phase  variations  across 
the  wavefront  remain  small,  i.e.,  <Ad>2>'^  less  than  1  radian.  Here,  the  weak  or  single¬ 
scatter  approximation  is  adequate  to  describe  the  resulting  amplitude  variations.  The  result 
is  that  a  set  of  receivers  at  some  distance  Z  measures  amplitude  variations  of  spatial  size 
about  equal  to  the  Fresnel  scale,  suitably  modified  by  the  propagation  geometry.  The  energy 
lost  by  the  main  wave,  which  goes  into  producing  the  amplitude  fluctuations,  is  small  and 
so  also  is  the  rms  variation  of  the  amplitude  fluctuations. 

As  the  wave  continues  to  propagate  and  A L  becomes  much  larger  than  the  outer 
scale  size  of  the  irregularities,  the  wave  is  “scattered”  many  times.  If  the  amplitude  of  the 
irregularities  is  large  enough,  this  results  in  an  rms  variation  of  phase  across  the  wavefront 
which  can  be  much  greater  than  unity.  Here,  then,  it  is  the  largest  scale  size  irregularities 
which,  after  many  scatters,  produce  the  largest  part  of  the  phase  variations.  When  <Ad>2>'* 
becomes  large  enough,  the  ionosphere  begins  to  act  like  an  optical  lens  and  produces  focus¬ 
ing.  Receivers  at  distance  Z  then  measure  amplitude  variations  which  are  smaller  in  spatial 
extent  than  the  Fresnel  scale.  These  small  scale  variations  are,  in  fact,  being  produced  by 
the  large  scale  irregularities  in  the  medium.  The  energy  lost  by  the  main  wave  in  producing 
these  amplitude  variations  is  considerable,  and  the  rms  variation  of  the  amplitude  fluctua¬ 
tions  can  be  very  large. 

From  the  above  description,  we  see  that  the  scattering  one  encounters  for  transio- 
nospheric  propagation  differs  according  to  various  situations.  If  the  extent  of  the  propaga¬ 
tion  region  is  small  or  the  amplitude  of  the  irregularities  is  also  small,  then  the  weak-scatter 
theory  may  be  used  to  describe  the  phenomena.  The  result  is  that  weak  amplitude  varia¬ 
tions  develop  which  have  a  characteristic  size  of  about  the  Fresnel  scale  for  the  receiving 
location.  However,  if  the  propagation  distance  extends  over  large  distances,  multiple  scat¬ 
tering  of  the  wave  can  produce  large  variations  in  the  rms  phase.  Here,  it  is  the  largest  size 
irregularities  which  act  to  produce  small  scale  amplitude  variations  at  a  receiver.  It  is  this 
latter  instance  of  strong  or  multiple  scattering  that  we  are  mainly  interested  in  in  this 
report.  In  the  following  section,  we  will  describe  this  scattering  process  in  more  quantitative 
detail  by  introducing  the  concept  of  the  phase-changing  screen. 
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MULTIPLE  REFRACTIVE  SCATTER  MODEL 

Before  we  can  understand  the  multiple  refractive  scatter  theory,  we  will  need  to 
describe  the  phase-screen  modeling  technique.  This  is  done  most  easily  by  using  the  concept 
of  an  angular  spectrum  of  plane  waves  (Booker  et  al .,  1950).  Such  a  spectrum  can  be  used 
to  synthesize  an  arbitrary  function  of  spatial  coordinates  in  terms  of  component  plane 
waves,  just  as  the  Fourier  integral  can  be  used  to  analyze  a  suitable  temporal  function  in 
terms  of  its  component  frequencies.  We  will  use  the  angular-spectrum  method  to  determine 
the  diffraction  field  produced  by  a  sinusoidal  phase-changing  screen.  Although  this  is  an 
extreme  simplification  of  our  actual  problem,  it  does  allow  the  generation  of  results  which 
are  also  valid  for  the  more  general  situation  to  be  described  later. 

For  this  section,  we  will  follow  Ratcliffe  (1956)  and  consider  plane  waves  normally 
incident  on  a  one-dimensional  phase  screen.  This  will  be  generalized  to  a  two-dimensional 
screen  and  obliquely  incident  spherical  waves  when  we  describe  our  implementation  of  the 
model. 

To  introduce  the  concept  of  an  angular  spectrum  of  plane  waves,  we  consider  an 
infinite  two-dimensional  «creen,  the  X-  Y  plane  at  Z  =  0,  which  is  uniform  in  the  Y  direc¬ 
tion.  Let  a  system  of  uniform  plane  waves  of  frequency  / and  wave  number  k  travelling  in 
directions  making  angles  +  and  -6  with  the  Z  axis  make  up  an  angular  spectrum  of  plane 
waves.  We  assume  that  the  waves  are  polarized  in  the  X-Z  plane  and  that  the  resultant 
complex  amplitude  of  those  waves,  with  wave  normals  between  8  and  dd ,  is  A(8)d8.  Then 
the  corresponding  bundle  of  waves  is  given  by 

A(8)d8  eikfX  sin  cos 


where  a  time  factor  e"'tu'  has  been  omitted.  The  total  field  amplitude  at  a  point  X,Z  is  given 
by 


fWi2  AmC,0,-S)e'Uxs*ZC)d8 

-w/2 

where  C  =  cos  8  and  5  =  sin  8.  In  particular,  the  amplitude  of  the  X  component,  call  it 
w(X,Z),  is  given  by 

w(X,Z)  =  /  n  2  A(8)Celk<xs*ZC)dd 

-tt/2 


or  with  d8  -  dS/  C  and  A(8)  =  W/(S') 


w(X,Z)=  /  1T(S)  eik{XS*ZCldS  (1) 

i 

This  represents  a  wave  field  varying  over  X  and  Z  written  in  terms  of  an  angular 
spectrum,  IF(5),  of  plane  waves.  Over  the  screen  at  Z  =  0  we  have 


w(A)  =  /  |  1T(5)  eikxsdS 
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If  we  let  k  -  2n/k  and  measure  distance  in  units  of  wavelength  k,  this  becomes 
w(x)  =  /  |  W{S)  ei2nSxdS 

where  x-X 'Ik.  For  the  remainder  of  this  section,  lower  case  letters  represent  distances  in 
units  of  wavelength  unless  otherwise  stated. 

From  these  results,  we  see  that  W^S)  and  w(x)  are  Fourier  transform  pairs  if  we 
extend  the  integration  limits  to  ±°°.  We  will  do  this  since  values  of  1 1  >  1  correspond  to 
evanescent  waves  which  decay  within  several  wavelengths  of  the  screen.  Consequently,  at 
planes  separated  from  the  screen  by  more  than  a  few  wavelengths,  only  propagating  fields 
for  which  1 5 1  <  1  will  exist.  So  we  can  write 


w(x)  =  /  "  fF(S)  eQnSxdS 
-00 


(2) 


and 


IF(S)=  /“  w{x)e-anSxdx 

_oo 


(3) 


Likewise,  Eq.  1  becomes 


w(x,z)  =  /  "  MS)  ean{xS*sCt  dS 

oo 


(4) 


The  above  analysis  shows  that  an  angular  spectrum  1F(S)  will  produce  on  the 
screen  at  Z  =  0  a  field  component  along  X  given  by  w(x),  where  w(jr)  and  W'XS)  are  trans¬ 
form  pairs.  This  implies  that  a  field  w(x)  over  the  screen,  however  it  may  be  produced,  is 
equivalent  to  the  angular  spectrum  1F(5). 

As  a  simple  example,  consider  an  angular  spectrum  consisting  of  just  two  plane 
waves  making  angles  ±0O  with  the  Z  axis  and  having  complex  amplitudes  1F(S0)  =  (a/ 2) 
exp  (id>)  and  kF(-50)  =  (a/2)  exp  (-id>),  where  S0  =  sin  60.  Then,  from  Eq.  4,  the  resulting 
field  amplitude  at  a  point  x,z  is 

w(x,z)=  lF(So)  ear(  s°x+c°z)  +  lF(-50)  e'2,H'V+Col) 

=  a  e,27rC°*  cos  (27rS0jr  +  d>) 

This  represents  a  waveform  travelling  along  the  +Z  direction  with  an  amplitude 
which  varies  sinusoidally  along  the  wavefront  with  wavelength  D.  The  spatial  period 
d-Djk  of  the  sinusoidal  variation  is  given  by  2nS0  -2tr  ox  d-  1  /  S0.  In  Fig.  2.  we  show 
the  resulting  wavefront  and  the  corresponding  angular  spectrum  (Ratcliffe,  1956). 

The  velocity  in  the  Z  direction  of  this  striated  wavefront  is  given  by  v0  =  //  C0, 
where  we  recall  that  we  are  assuming  sinusoidal  time  development  and  that  distances  are 
measured  in  units  of  k.  If  we  let  v  =  / be  the  velocity  of  a  uniform  plane  wave  and  write 
C0  =  (I  -  \j d2)'/\  then 
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Figure  2.  Illustrating  the  angular  spectrum  and 
resulting  waveform  for  a  two-component  system 
(after  Ratcliffe,  1956). 

v 

V°  "  (1  -  1  Id1)* 

The  velocity  of  such  a  wave  is  greater  than  that  of  the  uniform  plane  wave  and  depends  on 
the  spatial  period.  For  large  d  (small  S0),  the  velocity  differs  only  slightly.  However,  for 
small  d  (large  S0),  the  velocity  of  the  striated  wave  can  differ  significantly  from  the  uniform 
wave. 

On  the  plane  Z  =  0,  we  have 

w(x )  =  a  cos  (2tt50x  +  d>) 

This  field  could  have  been  produced  by  allowing  a  uniform  plane  wave  travelling 
along  0 Z  to  pass  through  a  thin  screen  which  imparts  this  amplitude  on  it.  The 
resulting  disturbance  would  be  the  same  as  that  produced  by  our  two  plane  waves 
at  0O  and  -0O.  They  thus  represent  the  angular  spectrum  corresponding  to  diffraction  by 
such  a  screen. 

The  power  of  this  representation  lies  in  the  fact  that  any  function,  regular  or 
random,  can  be  synthesized  from  an  angular  spectrum  of  plane  waves.  Thus  the  technique, 
described  here  for  a  regular  function,  can  be  used  for  a  randomly  varying  screen,  and  most 
of  the  results  we  derive  here  are  equally  true  for  that  more  general  case. 

Consider  now  a  thin  diffracting  screen  situated  in  the  Z  =  0  plane  which  alters  the 
phase  of  the  incident  wave  without  altering  its  amplitude.  This  approximation  is  valid  for 
F-region  scintillation  modeling  when  the  frequency  of  the  transmitted  wave  is  much  larger 
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than  the  ion-electron  collision  frequency.  We  suppose  that  as  the  wave  emerges  from  such 
a  screen,  the  phase  varies  along  the  X  direction  according  to 


<t>(x)  =  <Ad>2>^  cos  (27 rx/d) 

where  <Ad>2>'^  is  the  rms  value  of  the  phase  deviation  imparted  by  the  screen.  When  an 
“equivalent”  screen  is  used  to  replace  an  extended  propagation  region,  such  as  the  F-region 
of  the  ionosphere,  <Ad>2>'^  represents  the  rms  value  of  the  phase  fluctuations  measured  at 
a  receiver  after  the  wave  has  traversed  the  medium. 

For  a  normally  incident  uniform  plane  wave  of  unit  amplitude,  the  field  at  Z  =  0  is 
u’(x)  =  et'<A^2>^cos(2'rjf/^]  (5) 


In  order  to  determine  the  angular  spectrum  corresponding  to  Eq.  5,  we  use  the  familiar 
expansion,  with  r)  -  2nx/d  (Gradshteyn,  Ryzhik,  1965): 

e[ »< A02>’^cos( t,)]  =  (/nyn« Ad>2>'^)  timr' 

n =-» 

where  J  is  the  cylindrical  Bessel  function  of  the  first  kind. 

Thus 


w(x)  =  S°°  ( in)Jn[<A(t>2>,/! ]  e""> 

n=-°° 

We  showed  above  that  w(x)  and  ^(5)  are  Fourier  transform  pairs,  and  so  we  can 
write  Eq.  2  as 

w(x)  =  l°°  mSj  eanS"x 

n-  °o 

Comparing  these  two  expressions  for  w(x),  we  see  that  the  spectrum  has  components 
at  values  51  =  0, 1  /  rf,  - 1  /  d,  2/  d,  -2j  d, . . . . ,  n/d,  -n/d. . . .  The  amplitude  of  the  side  compo¬ 
nents  is  inJn  (<A</>2>,/5)  and  changes  by  7t/2  between  successive  components.  However, 

Jn  f<A d>2>l/j)  is  small  for  n  >  <A d>2>'^  and  so  there  is  little  energy  in  those  components 
for  which  S  >  <Ad>2>'/V^. 

For  <A d>2>'^  <  1  radian,  the  angular  spectrum  for  the  sinusoidal  phase  screen 
extends  to  components  at  1  Id  and  - 1  / d.  The  amplitudes  of  the  side  components  are  in 
quadrature  with  the  main  component  at  S  =  0.  If  <A <t>2>,/!  >  1  radian,  the  so-called 
“deeply  modulated”  phase  screen,  the  angular  spectrum  extends  to  components  at  approx¬ 
imately  +  and  -  <A <^2>'^/cf.  Each  component  amplitude  increases  tt/2  in  phase  relative 
to  the  undeviated  component  at  S'  =  0.  In  particular,  the  largest  angle  component  at 
5~<Ad>2>'/2/^ differs  in  phase  with  respect  to  the  S  =  0  component  by  <A</>2>'^7t/2.  The 
angular  spectrum  here  depends  jointly  on  the  spatial  period  and  the  depth  of  the  modula¬ 
tion  <Ad>2>'^. 

We  have  established  the  relationship  between  the  field  distribution  at  the  diffracting 
screen  and  the  angular  spectrum.  Now  we  want  to  determine  the  field  distribution  at  planes 
separated  from  the  phase  screen.  In  particular,  at  a  plane  located  at  some  z-z'  where 
z»  1. 
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Again  we  consider  a  uniform,  unit  amplitude  plane  wave  normally  incident  on  a 
one-dimensional  phase-changing  screen  at  Z  =  0.  The  wave  passes  through  the  screen,  and 
the  field  amplitude  at  any  point  x,z  is  given  by  Eq.  4.  As  the  striated  wave  formed  by  spec¬ 
tral  components  in  direction  ±5  propagates  from  a  point  x,0  to  a  point  x.f,  its  phase,  /3, 
changes  by  2tt(Sx  +  Cf)  -  InSx  or,  with  C  =  (1  -  S2)'A-  and  S«1 


A/3  =  2jtz’  -  nS2z’ 


The  first  part  of  this  equation  represents  the  phase  change  of  the  direct  ( S  =  0)  component 
in  travelling  to  the  z-z"  plane.  The  second  part  represents  the  phase  gain  by  the  striated 
wave.  If  g(x)  represents  the  complex  field  amplitude  on  the  plane  at  z-z!  with  phase  refer¬ 
enced  to  the  point  x,f ,  then  g(x)  is  synthesized  from  the  same  angular  spectrum  as  that 
which  makes  up  w(x)  on  the  phase  screen  at  z  =  0.  However,  each  component  wave  is 
advanced  in  phase  by  ttz’S2.  This  is  entirely  equivalent  to  the  usual  Huygen-Fresnel 
method  of  determining  the  diffraction  field  via  equivalent  current  sources  and  Greens  func¬ 
tion  (Ratcliffe,  1956). 

As  an  example  of  the  above  analysis,  we  apply  it  to  the  sinusoidal  phase  screen.  Let 
the  field  at  the  screen  be  given  by  Eq.  5.  If  <A<p2>'A  <  1  radian,  we  can  expand  to  first 
order  to  obtain 

w(x)  =  1  +  i<A<(>2>'A  cos  (27 rx/d) 
or,  equivalently 


i<Ad>2>'A 

w(jc)  =  1  +  - - -  (ei27rSx  +  c-i2n Sx) 

where  we  have  let  S  =  1/d.  Thus,  the  spectrum  contains  components  at  S  -  0,  \/d,  and 
-l /d,  as  we  indicated  earlier. 

It  was  shown  earlier  that  the  striated  wave  formed  by  the  two  side  components 
travels  at  a  greater  velocity  than  the  plane  wave.  At  a  certain  distance,  the  striated  wave, 
which  is  in  quadrature  with  the  main  wave  at  the  screen,  will  advance  in  phase  with  respect 
to  the  main  wave  and  will  be  in  phase  or  antiphase  with  it.  At  this  distance,  the  two  waves 
will  combine  to  produce  amplitude  variations  in  the  total  field.  From  above,  this  occurs 
at  the  plane  z-z'  where  irS2z'  =  v/ 2  or  z'  =  d2/ 2.  The  field  distribution  at  this  plane  is 
given  by 

g(x)  -  I  -  <A 4>2>'a  cos  ( 2nx/d) 

Thus,  the  period  of  the  amplitude  variation  is  the  same  as  that  of  the  phase  variation  in  the 
screen. 

As  the  waves  continue  to  advance  in  the  Z  direction,  the  striated  wave  gains  jt/2 
radians  in  phase  on  the  main  wave  at  each  distance  increment  of  d2/  2.  When  they  are  in 
phase  or  antiphase,  they  combine  to  produce  amplitude  variations,  and  they  produce  phase 
variations  when  they  are  in  quadrature. 

These  results  show  that  when  the  phase  fluctuations  are  weak  (<A</»2>W  <  1), 
amplitude  fluctuations  first  develop  on  a  plane  at  distance  T  from  the  diffracting  screen, 
for  which  F  -  D.  Here  D  is  the  spatial  period  of  the  phase  variations  at  the  screen  and 
F  -  \J2KZ’  is  the  Fresnel  zone  size  on  the  screen  at  a  distance  Z'.  The  amplitude  fluctua¬ 
tions  have  the  same  spatial  size  as  the  phase  variations  in  the  screen. 
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Next  we  consider  a  sinusoidal  phase  screen  for  which  <Ad>2>^  »  1  radian.  Here, 
the  angular  spectrum  extends  to  components  for  which  |5|~n/rf,  where  n~<A<£2>^.  Each 
of  these  side  components  is  advanced  by  n/2  radians  in  phase  on  its  predecessor  and  by 
nn/2  on  the  main  wave  (5  =  0).  If  we  examine  the  field  at  a  series  of  planes  at  increasing 
distances  from  the  diffracting  screen,  the  striated  waves  formed  by  component  pairs  will 
“catch  up”  with  the  main  wave  at  certain  distances.  The  component  in  the  direction  nS  will 
gain  n/2  radians  in  phase  on  the  main  wave  when  it  has  travelled  a  distance  nzn2S 2  =  nn/2 
or  z  -  (2  nS2)  1 .  At  this  distance,  the  component  wave  will  combine  with  the  main  wave 
to  produce  an  amplitude  modulation  with  a  spatial  period  of  \/ nS  or  d/ n.  As  the  wave 
advances,  each  successive  component  wave  will  combine  with  the  main  wave  to  produce  a 
different  amplitude  modulation  wavelength.  The  higher  order  component  waves,  corres¬ 
ponding  to  smaller  spatial  components,  travel  fastest  and  are  the  first  to  produce  amplitude 
modulation  at  planes  closest  to  the  diffracting  screen.  Thus,  at  all  distances,  the  diffraction 
field  contains  spatial  periods  smaller  than  the  spatial  period  of  the  screen  itself. 

To  see  these  effects  more  quantitatively,  consider  the  largest  order  component  wave 
for  which  S~<A<f>2>/! / d.  This  component  wave  produces  its  effect  on  the  amplitude  fluc¬ 
tuations  on  a  plane  at  z'  determined  by  nz'(<Atj>2>'/!  /  d)2  =  <Ad>2>'^  rr/2  or  for  z’  given  by 
z’  -  d2/ 2<Ad>2>'^.  In  terms  of  the  Fresnel  zone  size  F ,  this  becomes 

F  -  D/<A<i>2>'/,  (6) 

We  then  have  amplitude  variations  developed  on  planes  for  which  the  Fresnel  zone  is 
smaller  than  the  spatial  period  of  the  phase  screen.  This  is  in  contrast  to  the  results  for 
<A02>l/j  <  1,  where  we  found  amplitude  variations  first  develop  in  a  plane  for  which 
F  =  D.  The  large  phase  fluctuations  here  are  causing  the  screen  to  act  as  an  optical  lens  to 
produce  focusing  in  a  plane  closer  to  the  screen. 

The  spatial  period,  call  it  /,  of  the  amplitude  variations  produced  on  the  plane  at  z’ 
is  d/<A<t>2>l/! .  Then  with  /  =  L/A,  we  have 

L  =  D/<  A<£2>K  (7) 

and  we  see  that  the  smallest  period  amplitude  fluctuations  are  less  than  the  spatial  period 
of  the  phase  screen.  Again,  this  is  to  be  contrasted  to  <A <f>2>'/2  <  1,  where  we  found  the 
spatial  period  of  the  amplitude  variation  is  D. 

Now  consider  a  sinusoidal  phase  screen  at  Z  =  0  and  a  reception  plane  at  some 
distance  Z  for  which  the  Fresnel  zone  size  is  F.  Let  <A d>2>'*  »  1  and  D  »  F.  For  these 
conditions,  we  are  concerned  with  the  lens-like  refractive  scattering  (Booker,  MajidiAhi, 
1981)  described  above.  Such  a  screen  will  produce  foci  in  a  plane  parallel  to  the  screen. 
These  foci  lie  in  the  reception  plane  at  Z  if,  from  Eq.  6 

D  =  F<Ad>2>'/4  (8) 

Thus,  to  focus  a  sinusoidal  phase  screen  of  wavelength  D  at  a  reception  plane  specified  by 
F  requires  that 

<Ad>2>  =  (D/F)4  (9) 
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In  such  a  focal  plane,  the  wavelength  D  can  be  determined,  since  it  represents  the  distance 
between  focal  spots.  The  focal  spot  size  is  given  by  Eq.  7.  Thus,  for  a  given  screen  wave¬ 
length,  the  larger  the  rms  plane  fluctuations,  the  closer  the  focal  plane  to  the  screen  and, 
also,  the  sharper  the  focal  spot  size. 

The  fine  structure  in  the  reception  plane  is  associated  with  the  arrival  at  the  plane 
of  angular  spectrum  components  which  are  approximately  co-phased  within  a  cone  angle 
of  ±\/L.  This,  in  turn,  is  controlled  by  the  distance  across  the  phase  screen  over  which  the 
field  has  approximately  uniform  phase,  a  distance  given  by  Eq.  7.  So  for  refractive  scatter¬ 
ing,  fine  structure  in  any  reception  plane  is  the  same  as  it  is  in  the  focal  plane.  However, 
in  planes  displaced  from  the  focal  plane,  the  amplitude  variation  no  longer  has  the  form  of 
an  array  of  focal  spots  of  width  L  and  separation  D.  Out  of  the  focal  plane,  interference  by 
other  component  waves  with  slightly  different  phase  causes  an  amplitude  variation  with 
fine  structure  L  superimposed  on  a  gross  modulation  of  wavelength  D  (Booker,  MajidiAhi, 
1981). 

In  summary,  we  see  that  one  can  describe  the  effects  of  a  sinusoidal  phase-changing 
screen  in  terms  of  the  sizes  of  the  amplitude  variation  structure  which  it  produces  on  dis¬ 
placed  planes.  When  <A</>2>'^  <  1,  amplitude  variations  are  produced  at  distant  planes  for 
which  the  Fresnel  zone  size  is  equal  to  the  wavelength  in  the  screen.  The  amplitude  varia¬ 
tions  have  a  characteristic  size  approximately  equal  to  the  wavelength  in  the  screen.  This  is 
commonly  referred  to  as  weak  or  single  scattering  and,  for  the  ionosphere,  is  generally  the 
situation  for  very  high  frequencies  or  very  quiescent  conditions. 

In  contrast  to  the  weak  scattering  situation,  we  have  the  strong  or  refractive  scatter¬ 
ing  condition  for  which  <Ad>2>'/:  »  I  radian.  For  strong  scattering,  higher  order  compo¬ 
nents  in  the  angular  spectrum  produce  fine  amplitude  variation  structure  «D).  The  fine 
structure  is  produced  on  planes  closer  to  the  phase  screen. 

These  results  now  need  to  be  generalized  to  the  more  complicated  irregular  phase 
screen.  As  we  indicated  earlier,  the  irregular  structure  of  the  ionosphere  is  characterized  by 
a  random  distribution  of  irregularities  covering  a  large  range  of  scale  sizes.  In  order  to  use 
the  phase-changing-screen  modeling  technique,  we  must  relate  the  observed  irregular  struc¬ 
ture  of  the  density  irregularities  to  the  properties  of  the  phase  screen.  To  do  this,  we  will 
relate  the  power  spectrum  of  the  phase  variations  produced  by  the  screen  to  the  power 
spectrum  of  the  irregularities. 

Consistent  with  the  above  discussion,  we  consider  a  wave  which  is  normally 
incident  on  a  plane  ionospheric  layer  of  thickness  AZ..  We  assume  the  background  electron 
density  is  a  function  of  the  vertical  coordinate  z  only,  where  x,y,z  now  have  the  usual  units 
of  distance.  Within  the  plasma  layer,  there  exist  electron  density  irregularities  describable 
by  a  three-dimensional  power  spectrum,  S»  (kx,ky.k:).  If  we  consider  only  that  part  of 
the  spectrum  for  which  the  irregularities  are  much  larger  than  the  wavelength,  we  may 
apply  geometric  optics.  Under  this  approximation,  the  phase  at  a  point  x,y  in  a  plane  at 
z,  z  >  AL,  is  given  by 

<t>(x.v)  =  k  /  n(x,  v,z)dz  (10) 

o 

where  n  is  the  index  of  refraction  within  the  layer.  The  variation  in  the  phase  observed  at 
the  point  x,y  due  to  variation  in  the  index  of  refraction  along  the  ray  is 
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(II) 


c  M 

A<f>  ( x,y )  -  k  J  An(x,y,z)dz 

o 

Note  that  in  Eq.  1 1,  we  are  calculating  the  phase  variation  due  to  spatial  variation  in  the 
index  of  refraction  along  a  deterministic  ray.  Under  very  strong  scattering  conditions,  there 
may  be  many  rays  reaching  the  point  x,y.  In  those  conditions,  Eq.  1 1  represents  the  phase 
variation  along  a  “representative”  ray. 

If  we  do  a  similar  calculation  for  a  ray  which  reaches  a  point  x  +  Ax,  y  +  Ay  and 
form  the  averaged  product  with  Eq.  1 1,  we  have 


<A<l>(x,y)A<t>(x  +  A  x,y  +  Ay)> 


k2  /  A L  f  AL  <An(x,y,z)An(x  +  A x,y  +  A y.z'^dzdz’ 
o  o 


(12) 


where  we  assume  An  is  a  random  function  of  position  and  brackets  represent  the  ensemble 
average,  i.e.,  average  over  all  possible  configurations  of  the  random  function  at  a  given 
point. 

The  index  of  refraction  in  the  ionosphere  is  actually  a  function  of  space  and  time; 
however,  in  Eq.  1 1  and  12,  we  have  ignored  this  temporal  dependence.  The  assumption  is 
that  during  the  time  that  a  wave  takes  to  traverse  the  distance  A L,  the  ionosphere  is 
approximately  constant,  except  for  bulk  drift  motion.  This  assumption  of  a  frozen  iono¬ 
sphere  simplifies  the  calculations. 

To  proceed,  we  require  a  model  for  the  variation  in  index  of  refraction  appearing  in 
Eq.  12.  For  HF  and  higher  frequencies,  we  will  ignore  the  propagation  effects  of  the  earth’s 
magnetic  field  and  the  damping  effects  of  electron-ion  and  electron-neutral  collisions.  The 
ionosphere  is  then  describable  by  the  cold  plasma  dispersion  relation 


n  =  (1  -  e2 N / e0/mw2)'^ 


(13) 


where  e  and  m  are  the  charge  and  mass  of  the  electron,  respectively,  e0  is  the  free-space 
permittivity,  and  w  =  27r/is  the  operational  angular  frequency.  N(z)  is  the  height-dependent 
electron  density. 

For  our  purposes  of  relating  the  spatial  properties  of  the  irregularities  to  those  of 
the  phase  spectrum,  we  will  assume  <u 2  »  e2N/t0  m  and  so,  to  first  order  in  AN,  we  have 
from  Eq.  13 


An 


e2A  N 
2  e0mou2 


Assuming  AN  is  a  zero  mean,  homogeneous,  random  function  of  position,  Eq.  12  becomes 


p^Ax.Ay)  o c  J  J  pN(Ax,Ay.z  -  z’)  dzdz ’ 
0  0 


(14) 


where  p^,  pN  are  the  autocorrelation  function  of  phase  and  electron  density,  respectively. 
In  Eq.  14,  the  homogeneity  of  AN  has  allowed  us  to  write  the  autocorrelation  as  a  function 
of  the  difference  in  coordinates  between  the  two  points  (Panchev,  1971). 
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The  integral  in  Eq.  14  is  simplified  by  a  change  to  the  variables  a-z-z', 

/3  =  (z  +  z’)l  2.  Then,  for  a  layer  thickness  much  larger  than  the  outer  scale  of  irregularities, 
Eq.  14  becomes 


P*  (Ax, Ay)  oc  /  pN(Ax,Ay,a)da  (15) 

_oo 

The  integration  limits  have  been  extended  to  ±«  in  Eq.  15  since  pN  -  0  for  a  greater  than 
the  vertical  outer  scale  size.  This  equation  relates  the  two-dimensional  spatial  autocorrela¬ 
tion  function  of  phase  to  that  of  the  irregularities  within  the  layer.  The  horizontal  power 
spectrum  is  obtained  from  Eq.  15  by  Fourier  transformation  in  horizontal  coordinates 
(Panchev,  1971;  Tatarskii,  1971). 

For  a  three-dimensional  autocorrelation  function  of  irregularities,  the  power  spec¬ 
trum  SN  is  given  by 


SN  (kx,ky,k,)  =  /  /  /  pN(Ax,Ay,Az)  tik  Ar  d(Ar) 


where  A r  -  (Ax.Ay.Az). 

A  two-dimensional  spectrum  in  the  plane  at  z  is  then  given  by 
Fn  (kx.ky,Az)  -  f° J  p N(Ax,Ay,Az)  ei(k*Ax'k>Av) d(Ax)d(Ay) 


From  these  two  expressions,  we  see  that 


SN  (kx,kv,kz)  =  r  Fn  (k x,ky,A:)  eik:A:d(Az) 

oo 


So,  if  we  perform  a  two-dimensional  Fourier  transform  of  Eq.  15  with  respect  to  Ax,  A.y, 
we  find 


<t>(kx,kv)ccSN(kx,ky.O)  (16) 

Thus,  the  two-dimensional  power  spectrum  of  phase  in  the  plane  at  z  is  simply 
related  to  the  three-dimensional  spectrum  of  irregularities.  In  particular,  for  a  power  law 
irregularity  spectrum,  the  phase  spectrum  is  also  of  power  law  form.  Integration  of  Eq.  16 
along  y  shows  that,  similarly,  the  one-dimensional  phase  spectrum  is  proportional  to  the 
two-dimensional  irregularity  spectrum. 

Uscinski  el  al.  (1981)  have  treated  a  one-dimensional  phase  screen  with  a  power 
law  spectrum  of  phase  fluctuations  with  spectral  index  4.  This  corresponds  to  a  two- 
dimensional  screen  with  spectral  index  5.  When  the  rms  phase  fluctuations  are  large 
compared  to  I  radian  and  the  outer  scale  of  irregularities,  La,  is  large  compared  with  the 
Fresnel  scale,  Uscinski  el  al.  find  that  the  fine  scale  structure  in  the  intensity  variation 
pattern  is  given  approximately  by 


/  = 


<A<t>2>,/! 


(17) 
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Comparison  with  Eq.  7  shows  that  in  going  from  a  one-dimensional  sinusoidal 
phase  screen  to  an  irregular  power  law  screen  with  spectral  index  4,  the  fine  structure,  or 
focal  scale,  is  determined  by  replacing  the  sinusoidal  wavelength  by  the  outer  scale  of  irreg¬ 
ularities  (Booker,  MajidiAhi,  1981). 

Booker  and  Tao  (1987)  have  shown  that  for  a  one-dimensional  phase  screen  with 
spectral  index  I  <p<3,  Eq.  17  must  be  modified  to 


(18) 


where  a  is  a  numerical  coefficient  which  is  a  function  of  p.  Similar  results  apply  to  a  two- 
dimensional  screen  with  a  spectral  index  p  +  1. 

In  Fig.  3,  several  examples  of  intensity  power  spectra  are  shown.  These  data  were 
derived  from  the  signal  of  the  ATS-6  beacon  satellite  at  Boulder,  Colorado,  in  January 
1975.  Simultaneous  data  at  three  frequencies  are  shown  in  the  figure.  At  the  two  higher 
frequencies,  the  S4  index,  which  is  the  normalized  variation  of  the  intensity  about  the 
mean,  is  less  than  1.  This  is  representative  of  weak-to-moderate  scattering.  At  40  Mhz, 
however,  the  S4  index  is  greater  than  I,  indicative  of  a  strong,  multiple-scatter  condition, 
as  described  above.  Note  that  the  spectra  in  this  case  extend  to  high  normalized  frequen¬ 
cies,  corresponding  to  the  fine  scale  structure  produced  under  these  conditions. 


Figure  3.  Illustrating  intensity  power  spectra  for  40, 

140,  and  360-MHz  signals  transmitted  from  the  beacon 
satellite  ATS-6.  The  140  and  360-MHz  spectra  are  Indi¬ 
cative  of  weak  scattering,  while  the  40-MHz  spectrum  is 
typical  of  strong  scattering.  (After  Umeki,  eta!.,  1977.) 
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Using  the  theory  outlined  above  for  the  deeply  modulated  phase  screen,  Booker 
and  MajidiAhi  (1981)  have  generated  the  intensity  spectra  shown  in  Fig.  4.  They  assumed 
a  power  law  phase  spectrum  with  a  spectral  index  p  -  3,  an  outer  scale  /  0  =  10  F,  and  an 
inner  scale  L,  =  10~2  F ,  where  Fis  the  Fresnel  scale.  Various  values  of  the  mean  square 
phase  fluctuation  were  used  in  Fig.  4  to  produce  the  multiple  plots. 

Referring  to  Fig.  4,  we  see  that  for  <Ad>2>  <  10  radians,  the  scatter  is  weak  and 
the  normalized  spectra  show  a  single  peak  at  approximately  the  Fresnel  scale.  As  <Ad>2>  is 
allowed  to  increase,  the  small  scale  structure  begins  to  emerge  and  the  spectra  extend  to 
higher  frequencies.  In  this  instance,  the  focal  scale,  /,  as  defined  above  and  indicated  by  the 
x  in  Fig.  4,  corresponds  approximately  to  the  high-frequency  roll-off  point  in  the  spectrum. 
Note  in  Fig.  4  that  the  lens  scale  is  also  indicated  on  the  plots.  Comparison  to  Fig.  3  shows 
that,  qualitatively,  the  predicted  spectra  are  similar  to  those  seen  under  strong  multiple- 
scatter  conditions. 

So  far,  we  have  indicated  that  the  rms  phase  fluctuation  must  be  larger  than 
1  radian  for  multiple  refractive  scattering  to  dominate  weak  diffractive  scattering.  However, 
in  Fig.  4,  weak  scatter  dominates  even  for  rms  phase  fluctuations  greater  than  3  radians. 

We  want  now  to  determine  approximately  how  large  the  phase  fluctuation  must  be  for  the 
small  scale  structure  given  by  Eq.  18  to  develop. 
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Figure  4.  Intensity  power  spectra  generated  by  using  the  multiple 
refractive  scatter  model.  Calculations  utilized  a  p  =  3  power  law 
phase  spectrum  with  an  outer  scale  L0  =  10F  and  an  inner  scale 
L ,  =  10”2F,  where  F  is  the  Fresnel  scale.  Curves  shown  are  para¬ 
meterized  by  the  mean  square  phase  fluctuations.  (After  Booker, 
MajidiAhi,  1981.) 
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Consider,  then,  a  random  plasma  slab  containing  irregularities  of  electron  density 
describable  by  a  power  law  spectrum  between  some  outer  scale  L0  and  inner  scale  L{  .  Also 
assume  that  a  reception  plane  parallel  to  the  slab  lies  at  some  distance  z.  For  most  iono¬ 
spheric  scintillation  phenomena,  L0  »  F ,  where  F  is  the  Fresnel  scale  for  the  plane  at  z. 

In  fact,  this  is  a  requirement  for  the  validity  of  the  multiple  refractive  scatter  theory. 

We  wish  to  model  this  situation  with  an  equivalent  phase-changing  screen  which 
will  replace  the  actual  scattering  medium.  Under  conditions  of  multiple  refractive  scatter¬ 
ing,  there  are  multiple  rays  joining  a  transmitter  to  a  receiver  in  the  reception  plane.  This  is 
illustrated  schematically  in  Fig.  5a,  with  the  equivalent  phase  screen  depicted  in  Fig.  5b. 
Associated  with  each  of  these  rays  is  a  Fresnel  zone  of  scale  F.  These  zones  are  distributed 
across  the  phase  screen  with  a  scale  L,  which  we  have  called  the  lens  scale,  given  here  by 

L=F2H  (19) 

where  /  is  the  focal  scale. 

In  order  for  there  to  be  many  Fresnel  zones  across  the  screen,  we  must  have 

L  »  F 
or ,  from  Eq.  19 

l«F  (20) 

From  Eq.  18,  then,  for  a  one-dimensional  screen  with  spectral  index  1  <p  <  3  or  a 
two-dimensional  screen  with  spectral  index  2<p<4,  multiple  refractive  scattering  domi¬ 
nates  diffractive  scattering  when 

<A  d>2>»i^yl  (21) 

Since  L0  »  Fand  p  >  1,  <A$2>  must  be  very  large  for  the  simple  theory  outlined 
above  to  hold,  or  equivalently,  the  wavefront  must  be  scattered  many  times. 

Now  the  focal  scale,  /,  is  associated  with  the  arrival  at  the  reception  of  rays  spread 
over  a  cone  of  half-angle  KI2ttI.  This  angle,  in  most  cases,  is  small  compared  with  1  radian. 
Under  severe  scattering  conditions,  however,  the  accumulated  scattering  angle  may 
approach  1  radian.  Under  such  conditions,  the  correlation  scale  in  the  reception  plane 
approaches  k/2n  and  the  ionosphere  is  acting  as  a  diffuse  scatterer.  This  is  the  minimum 
correlation  scale  observed  (Booker,  Tao,  1987). 

To  summarize,  the  correlation  structure  of  the  intensity  variation  varies  according 
to  Eq.  18.  When  /,  calculated  according  to  Eq.  18,  is  greater  than  the  Fresnel  scale  F,  the 
scattering  is  weak  and  intensity  variations  of  scale,  F,  are  observed  at  the  reception  plane. 
When  /  is  less  than  Fbut  greater  than  A/2?r,  multiple  refractive  scattering  is  dominating 
single  diffractive  scattering  and  the  intensity  correlation  scale  is  given  by  /.  When  /  <  \/2tt, 
the  ionosphere  has  become  a  diffuse  scatterer,  with  incoming  radiation  spread  over  a  cone 
half-angle  of  approximately  1  radian.  In  this  case,  the  correlation  scale  /  is  about  A/2t r. 

Using  the  multiple  refractive  scale  theory,  it  is  possible  to  determine  several  other 
properties  of  the  scattered  wave  which  are  useful  in  many  applications.  These  include  the 
correlation  bandwidth  (Booker,  Tao,  1987),  the  larger  scale  intensity  variation  at  the  lens 
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scale,  L,  which  Booker  called  “twinkling”  (Booker,  Tao,  1987),  and  the  half-angle  of  the 
cone  over  which  radiation  arrives  at  the  reception  plane  under  various  scattering  condi¬ 
tions.  For  our  purpose  in  this  report,  we  will  not  deal  directly  with  these  quantities.  The 
interested  reader  is  referred  to  the  Reference  section  for  more  information  (Booker,  1979; 
Booker,  Tao,  1987). 


Incident 

wavefront 


Scattering 

layer 


Figure  5.  Schematic  of  the  phase-screen  modeling 
method.  In  (b),  the  relationship  between  the  lens 
scale,  L,  and  the  Fresnel  scale,  F,  is  shown  for  a 
strong  scattering  condition. 
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IMPLEMENTATION  OF  THE  MULTIPLE 
REFRACTIVE  SCATTER  MODEL 

A  typical  operational  scenario  is  depicted  in  Fig.  6.  In  this  figure,  a  satellite  in  earth 
orbit  transmits  a  signal  at  low  VHF  or  HF  which  is  received  at  one  or  more  earth-based 
locations.  We  wish  to  determine  the  phase  correlation  of  the  signal  at  points  separated  in 
the  reception  plane  by  a  lag  Ar0. 

As  we  showed  earlier,  the  two-dimension  phase  power  spectrum  <t>(Jc)  in  a  plane 
parallel  to  the  equivalent  phase  screen  is  simply  related  to  the  three-dimension  power  spec¬ 
trum  of  electron  density  irregularities  SN(lc,k,)  by  the  relation 

m  oc  SN(Jc,k:  =  0) 

where,  referring  to  Fig.  6, k  =  kx  x  +  kyy.  In  our  current  version  of  the  model,  we  assume 
that  the  electron  density  irregularities  are  three-dimensionally  isotropic  with  a  power  law 
spectrum  for  irregularities  smaller  than  the  outer  scale,  L0  .  We  take  the  inner  scale  to  be 
zero. 

We  assume  initially  that  the  spectrum  is  adequately  described  by  a  single  spectral 
index,  p.  However,  as  we  indicated  earlier,  there  is  some  evidence  that  a  two-spectral-index 
model  may  be  necessary  in  some  situations.  Following  Booker  and  Tao  ( 1 987),  who  assumed 
a  one-dimensional  spectrum  with  spectral  index  of  2.4,  we  take  p  to  have  the  value  3.4. 

Thus,  we  write  the  phase  power  spectrum  in  a  plane  parallel  to  the  phase  screen  as 

<ttkx,ky)  =  <A</>2>*  S(kx,ky)  (22) 

where 


S(k*’ky)  °  [1  +  (k2  +  ky2)L02)Pl2 

<A<*>2>  is  the  mean  square  fluctuation  of  phase  and  L0  is  the  outer  scale  at  the  height  of 
the  phase  screen. 

Since 

— /  <t>(k)dk  =  <A d>2> 

(2  tt)z  °° 
we  must  have 


1  r  +«> 

(27t)2 


S(k)  dk  =  1 


So,  normalizing  S(k),  we  find 


S(kx,ky) 


4rrL02  T(p/2)  1 

T(pt  2-1)  [1  +  k2L02y12 


where  T  is  the  gamma  function  and  k 2  =  k2  +  k2. 


(23) 
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Figure  6.  Schematic  of  a  typical  operational  scenario,  indicated  is  the 
geometrical  relationship  between  the  phase  screen  and  the  transmitter 
and  receiver.  P,  and  P2  are  group  paths  along  the  indicated  portions 
of  the  ray,  6  is  the  zenith  angle  of  the  ray  at  the  reception  plane,  V  and 
tf>  are  the  zentih  and  azimuthal  angles  at  the  phase  screen,  respectively. 
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The  two-dimensional  phase  autocorrelation  function  in  a  plane  parallel  to  the 
screen  is  given  by  the  Fourier  transformation  of  the  power  spectrum.  So 

<A <P2>  Plt>  (AT)  =  — l—,  f°°  /  “  e‘*'7(t>0c)d]c  (24) 

( Z7T )z  -00 

when  AT  =  Axx  +  Ayy.  As  we  indicated  earlier,  we  are  assuming  the  fluctuation  of  electron 
density  is  describable  as  a  homogeneous  random  function.  In  that  case,  the  phase  auto¬ 
correlation  is  a  function  of  the  difference  in  coordinates  only. 

Combining  Eq.  22,  23,  and  24,  we  find 
2  C1 

P<t>(Ar)  =  —  (Arl1L0)  2  K (Ar/ L0)  (25) 


where  Ar 2  =  A*2  +  Ay2  and  Kv  is  the  modified  Bessel  function  of  the  second  kind. 

Referring  to  Fig.  6,  we  see  that,  in  general,  the  reception  plane  is  not  parallel  to  the 
phase  screen.  Therefore,  we  must  express  the  lag  A70  in  the  (x,y,z)  coordinate  system  of  the 
phase  screen.  Let  A70  have  components  (Ax,Ay,Az)  in  the  (x,y,z)  coordinate  system.  Then, 
assuming  the  rays  are  approxi  "Uely  straight  in  the  region  near  the  phase  screen,  we  have 

Axs  *»  Ax  -  tan  cos  d>  A z 

Ayx  ^  Ay  -  tan  i/r  sin  4>  A z 

where  (Axs,Ays)  are  the  components  of  the  equivalent  lag  in  the  phase  screen.  Here,  1 (1  is 
the  zenith  angle  of  the  ray  connecting  the  transmitter  to  the  receiver,  measured  at  the  phase 
screen,  and  d>  is  the  azimuth. 

Finally,  the  effective  lag  must  be  modified  to  include  the  spreading  effect  due  to  the 
finite  distance  to  the  point  source.  Referring  to  Fig.  6,  let  Pt  be  the  group  path  along  the 
ray  connecting  the  transmitter  to  the  penetration  point  in  the  screen  and  let  P2  be  that 
connecting  the  screen  to  the  receiver.  Then,  the  effective  lag  in  the  screen  corresponding  to 
A70  is  given  by 


P\ 

Pi*  Pi 


I  AT  s  | 


where 


(26) 


I  ArJ  =(A*t2  +  A  y2)* 

Perhaps  the  most  important  parameter  for  determining  the  nature  of  the  scattering 
for  transionospheric  propagation  is  the  mean  square  phase  fluctuation,  <A <f>2>.  As  we 
pointed  out  earlier,  phase  fluctuations  build  up  as  a  wave  propagates  through  an  extended 
medium.  A  single  receiver  also  measures  increased  phase  fluctuations  as  the  time  interval 
over  which  phase  measurements  are  made  is  increased.  This  is  due  to  the  effective  drift 
motion  of  the  ionosphere,  which  brings  larger  scale  irregularities  into  play  as  the  measure¬ 
ment  time  is  increased.  Phase  fluctuations  peak  when  sufficient  time  has  allowed  outer- 
scale  irregularities  to  contribute  to  these  fluctuations.  While  diffractive  scattering  by 
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smaller  scale  irregularities  (~  Fresnel  scale)  contributes  to  the  total  phase  fluctuations, 
most  of  the  content  is  due  to  refractive  scattering  by  outer-scale  irregularities.  This  is  true 
even  when  intensity  fluctuations  are  small. 

Consider  a  wave  incident  on  an  isotropic  electron  density  irregularity  of  scale  size 
L0  .  According  to  geometrical  optics,  the  mean  square  phase  fluctuation  due  to  such  an 
irregularity  is 

<A<(>2>0  =  k2La 2  <An2> 

where  <An2>  is  the  mean  square  fluctuation  of  index  of  refraction  associated  with  the 
irregularity.  Referring  to  Fig.  6,  let  As  be  a  segment  of  the  ray  connecting  PT  to  Pr  where, 
strictly  speaking,  As  »  LQ.  Then  there  are  approximately  As /  L0  irregularities  along  that 
ray  path  segment,  and  the  total  mean  square  phase  fluctuations  for  that  segment  is 

<A(f>2>  =  —  <Ad>2>0 

1-0 

The  total  mean  square  phase  fluctuation  along  the  ray  from  PT  to  Pr  is  then 

<A <t>2>  =  k2  /  Pj  La(z)  <A n2>  dJ  (27) 

R r 

where  integration  is  along  the  deterministic  ray  connecting  PT  to  Pr.  The  outer  scale  is  a 
function  of  height  and,  following  Booker  and  Tao  (1987),  has  the  profile  shown  in  Fig.  7. 
This  profile  approximates  the  variation  of  neutral  particle  scale  height  up  to  approximately 
400  km  (Booker,  Tao,  1987). 

In  order  to  calculate  <A </>2>  according  to  Eq.  27,  we  must  first  determine  the  ray 
through  the  background  ionosphere  which  connects  PT  to  Pr.  As  before,  we  assume  the 
background  ionosphere  is  describable  by  the  cold  plasma  dispersion  relation 


f  2 


N(r) 

AL 


(28) 


In  this  form,/p  is  the  peak  plasma  frequency  and  Nm  is  the  peak  electron  density.  In  this 
initial  version  of  the  model,  we  will  follow  Booker  and  Tao  (1987)  and  assume  that  N/  Nm 
is  a  function  of  height  only.  The  F-region  profile  we  will  use  is  shown  in  Fig.  8.  This  profile 
is  describable  by  the  functional  form  (Booker,  Tao,  1987) 


N 

~nZ 


-  exp 


-  Bln 


1  +  exp  [(z  -  zm)/C] 
1  +  exp  D 


(29) 


where  A  -  40  log  e  (km),  B  -  3/log  e,  C  =  100  km,  and  D  -  log  (C/ AB  -  E).  We  take  zm, 
the  height  of  the  peak  electron  density,  to  be  400  km. 


The  ray  connecting  PT  to  Pr  through  this  ionosphere  is  found  by  numerical  integra¬ 
tion  of  the  spherical  coordinate  ray  equation  (Kelso,  1964): 


X  =  re 


dr 

( n2r 4  -  r2r 2  sin  28)'/t 


(30) 
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Figure  7.  Variation  of  outer  scale  with  height 
used  in  calculation  of  the  mean  square  phase 
fluctuations. 


Figure  8.  Normalized  F-region  electron 
density  profile  used  in  the  current  version 
of  the  model. 
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Here  x  is  the  angle  subtended  at  the  center  of  the  earth  by  the  great  circle  path  from  the 
receiver  to  the  subtransmitter  point.  The  radial  coordinate  is  given  by  r  =  re  +  z,  where  re  is 
the  earth’s  radius  and  z  is  height  above  the  earth.  The  angle  0  is  the  angle  of  arrival  at  the 
earth-based  receiver  (see  Fig.  6).  The  ray  is  found  by  varying  the  angle  0  until  the  launched 
ray  meets  the  transmitter. 

It  remains,  in  Eq.  27,  to  determine  the  mean  square  fluctuation  of  index  of  refrac¬ 
tion.  Taking  differentials  in  Eq.  28,  we  have 


2nAn  =  - 


II 

f 


AN 


Squaring,  we  have 


1  fP2  &N2 
(2  n)2  f1  Nm2 


or,  taking  the  ensemble  average 


<A<f»2>  = 


(1  -  n2)2 

(2  n)2 


< 


(AN)2 

N„ 


> 


(31) 


The  procedure  for  determining  the  mean  square  phase  fluctuation  is  to  first  deter¬ 
mine  the  ray  path  through  the  background  ionosphere  which  connects  the  transmitter  to 
the  receiver.  This  is  done  by  iteration  on  6  in  Eq.  30.  Once  this  ray  is  found,  we  integrate 
Eq.  27  along  that  ray. 

Another  parameter  of  interest  is  the  mean  square  fluctuation  of  the  elevation  angle 
0  and  azimuthal  angle  r?  at  the  reception  plane.  For  a  one-dimensional  ionosphere  which 
we  are  currently  using,  the  rays  are  confined  to  vertical  planes  containing  the  transmitter 
and  receiver.  Thus,  we  choose  that  plane  as  the  azimuthal  reference  plane  and  so  rj  =  0  in 
the  following  development.  However,  since  our  phase  screen  is  two-dimensional,  we  can 
calculate  the  fluctuation  in  azimuth  around  zero. 

Let  3  be  the  vector  separating  two  receive  antennas  located  in  the  reception  plane 
(see  Fig.  6).  Let  d ,  be  the  component  of  3  in  the  vertical  plane  containing  the  receiver  and 
transmitter.  We  assume  d{  <  L0  where  L0  is  the  outer  scale  of  irregularities  at  the  phase 
screen. 

Let  <f>i jk  and  <t>i/k  be  the  phase  path  length  along  the  rays  connecting  the  trans¬ 
mitter  to  the  respective  receive  antennas.  Then  the  elevation  angle  0  in  the  reception  plane 
is  given  approximately  by 


sin  0  ** 


<t>  i  -  d>2 
kdx 


For  small  deviations  in  0,  we  have 


A<f>  i  -  A<t>j 

A(sin  0)  <**  A0  =  - 

kd  | 
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and  so  the  mean  square  fluctuation  in  elevation  angle  of  arrival  is  approximately 


<A02>  «  < 


(Ad»,  -  A<*>2)2  ^ 
k2d,2 


The  phase  structure  function  appearing  in  the  numerator  can  be  expanded  and,  assuming 
that  the  phase  autocorrelation  function  exists  and  that  the  statistics  of  the  medium  are 
homogeneous,  we  find 


<A02>  - 


2<Ad>2> 

k2dt2 


tl  -P^,)] 


(32) 


where  (d^  is  the  phase  autocorrelation  function.  As  in  our  earlier  discussion  concerning 
the  phase  autocorrelation,  the  effective  value  of  dt  ,  adjusted  for  a  tilted  reception  plane 
and  a  finite  source  distance,  must  be  used  in  Eq.  32.  If  d2  is  the  component  of  3  in  the 
direction  of  increasing  azimuth,  a  similar  equation  holds  for  the  fluctuations  <Atj2>  in 
azimuth,  with  d2  replacing  d\ . 

The  expression  for  the  mean  square  phase  fluctuation  appearing  in  Eq.  27  refers  to 
the  total  value  of  the  phase  fluctuation.  That  is,  this  is  the  value  one  would  measure  if  the 
signals  were  sampled  over  a  long  enough  period  of  time  for  outer-scale  irregularities  to 
contribute  to  the  fluctuations.  If  the  sampling  period  is  less  than  this,  the  outer-scale  irregu¬ 
larities  do  not  contribute  to  the  measured  phase  fluctuations  and  <A d>2>  is  reduced 
accordingly. 

For  the  measurement  of  the  signal  from  a  geostationary  satellite,  for  example,  the 
largest  scale  size  which  contributes  to  the  phase  fluctuation  is  determined  by  the  sampling 
period  and  the  drift  velocity  of  the  irregularities  at  the  height  of  the  phase  screen.  For  a  low 
orbiting  satellite,  the  drift  velocity  of  the  irregularities  is  normally  much  less  than  the  scan 
velocity  of  the  line-of-sight  between  receiver  and  transmitter.  In  this  case,  it  is  the  scan 
velocity  which  provides  an  effective  drift  velocity  at  the  phase  screen.  In  either  case,  if  the 
product  of  the  sampling  period  and  effective  drift  velocity  is  smaller  than  the  outer  scale, 
the  mean  square  fluctuation  of  phase  is  reduced. 

However,  it  is  also  true  that  it  is  not  always  necessary  to  sample  over  a  long  period 
to  bring  all  irregularity  scales  into  play.  Under  conditions  of  very  strong  scattering,  spread- 
F  conditions  for  example,  the  lens  scale  extends  over  many  Fresnel  scales.  A  reasonable 
sample  of  various  realizations  of  the  irregularity  distribution  can  then  be  made  over  a 
sampling  period  of  only  several  seconds  (Booker,  Tao,  1987). 

Thus,  it  is  necessary  to  adjust  the  mean  square  fluctuations  according  to  the  lens 
scale,  calculated  according  to  Eq.  18  and  19.  Booker  and  Tao  (1987)  have  shown  that  the 
effective  mean  square  phase  fluctuations  are  given  by 


<A  *2>f//  = 


[F/L0]p <Ad>2>,  L  <  F 
[L/L0\p-'  <Ad>2>,  F <  L<  La 
<A  <t>2>,  L<La 


(33) 


where  <A<£2>  is  given  by  Eq.  27.  It  is  this  value  of  the  mean  square  phase  fluctuation 
which  is  to  be  used  in  Eq.  32. 
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This  completes  the  description  of  the  model  currently  being  used  for  the  analysis 
of  angle-of-arrival  data.  Several  aspects  need  to  be  improved  for  the  accurate  testing  of  the 
method.  First,  a  more  realistic  ionospheric  model  must  be  used  in  conjunction  with  the 
irregularity  model.  It  is  expected  that  during  the  coming  year,  the  Ionospheric  Conductivity 
and  Electron  Density  (ICED)  model  (Tascione  et  al.,  1987)  will  be  incorporated  along  with 
a  full  three-dimensional  ray  trace  capability.  This  should  greatly  improve  our  ability  to 
accurately  determine  the  ray  connecting  the  transmitter  to  the  receiver. 

Further  work  is  also  needed  to  describe  the  properties  of  the  irregularities.  The 
spatial  and  spectral  properties  of  irregularities  of  scale  size  of  several  kilometers  is  fairly 
well  documented  and  some  modeling  does  exist  (Secan,  Fremouw,  1983).  However,  at  the 
lower  frequencies  of  interest  here,  a  major  role  is  played  by  large  scale  irregularities.  Very 
little  is  known  about  the  properties  of  these  irregularities,  since  their  effect  has  been  routinely 
removed  before  analysis  in  previous  studies.  As  we  indicated  in  an  earlier  section  of  this 
report,  evidence  indicates  that  the  spatial  gradient  of  electron  density,  AN,  is  approxi¬ 
mately  proportional  to  the  ambient  density,  so  that  AN/N  is  approximately  constant.  This 
assumption  is  made  in  our  current  model  and  is  used  in  our  calculation  of  the  rms  phase 
fluctuations.  As  crude  as  this  assumption  may  seem,  it  reflects  current  knowledge  in  the 
behavior  of  the  irregularity  amplitudes. 

The  three-dimensional  structure  of  the  large  scale  irregularities  is  also  an  area  in 
which  little  is  known.  Our  assumption  of  isotropic  large  scale  irregularities  represents  a 
theoretical  “best  guess”  based  on  the  mechanism  for  their  production. 

What  will  be  required  to  answer  these  questions  is  a  coordinated  experimental 
effort  designed  exclusively  to  investigate  the  nature  of  large  scale  irregularities.  Some  effort 
along  these  lines  has  already  begun,*  but  a  complete  global  description  will  require  a  much 
greater  effort. 


SUMMARY 

In  this  report,  we  have  indicated  methods  for  quantitatively  determining  the 
effects  of  strong  scatter  on  locating  systems.  The  theory  outlined  is  relatively  simple  and, 
in  conjunction  with  an  accurate  orbital  prediction  model,  should  provide  a  useful  tool  for 
analyzing  angle-of-arrival  data. 

Two  things  are  required  to  achieve  improvement  in  the  accuracy  of  the  method. 
First,  a  more  realistic  description  of  the  background  ambient  ionosphere  is  needed  for 
determining  the  ray  connecting  the  transmitter  to  the  receiver.  The  accurate  determination 
of  this  ray  will  provide  an  improved  estimate  of  the  actual  angular  location  of  the  transmit¬ 
ter,  i.e.,  the  magnitude  of  the  refractive  error  correction.  Also,  as  we  have  indicated,  the 
determination  of  the  fluctuation  in  this  angular  location  caused  by  scattering  from  irregu¬ 
larities  requires  the  integration  of  the  irregularity  spectrum  along  this  ray.  Thus,  a  more 
accurate  ray  determination  will  result  in  more  accurate  estimates  of  the  error  in  location 
due  to  scattering. 


*S.  Basu,  Emanuel  College,  Boston,  private  communication.  1988. 
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The  requirement  for  a  more  accurate  deterministic  ray  determination  should  be 
satisfied  by  the  use  of  the  ICED  electron  density  model  along  with  a  full  three-dimensional 
ray-tracing  capability.  This  improvement,  expected  shortly,  should  greatly  increase  the 
accuracy  of  both  aspects  of  the  problem. 

The  second  requirement  for  improved  accuracy  is  a  more  realistic  model  of  the 
irregularity  structure.  This  information  is  required  for  the  accurate  determination  of  the 
mean  square  phase  fluctuations  which,  as  shown  in  this  report,  is  basic  for  estimating  the 
magnitude  of  the  effect  of  the  scattering.  Work  along  these  lines  this  year  is  expected  to 
increase  the  accuracy  of  the  prediction  capability  of  this  model. 
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